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Abstract 

The Hanbury-Brown-Twiss (HBT) radii have been calculated from the two particle 
correlation functions with virtual photons produced in the collisions of two nuclei 
at ultra-relativistic energies. We show that the variation of the HBT radii with the 
invariant mass of the virtual photon can be used to characterize and distinguish 
the hadronic as well as the partonic phase that might have produced initially in the 
collisions. It has been illustrated that the non-monotonic variation of the HBT radii 
with invariant mass provides an access to the development of collective flow in the 
system. 
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The aim of the nuclear collisions at ultra-relativistic energies is to create and study 
- a state of matter, where quarks and gluons are dislocated from individual hadrons 
and make an excursion over a nuclear volume - such a phase of matter is called quark 
gluon plasma (QGP). Several probes - both electromagnetic (EM) and hadronic have 
been proposed for the diagnostics of QGP. The electromagnetically interacting probes 
(real and virtual photons) has the advantage over the hadronic probes because of 
their nature of interaction. The EM probes has mean free path much larger than the 
size of the system as a consequence they leave the system without re-scattering and 
hence can transmit the source information very efficiently [1]. However, photons and 
lepton pairs can be produced from both the partonic as well as hadronic phases [2]. 
As a consequence the disentanglement of the contribution for the QGP still remains 
a big challenge. 

In case of EM probes- dilepton has the advantage over the real photons. The photons 
with low transverse momentum (fc^-) from the hadronic phase may receive large 
transverse kick due to radial flow and consequently appear as high kr photons. 
These photons can mingle with the the contributions from the high temperature QGP 
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phase, making the detection of photons from QGP difficult. However, for dileptons 
there are two kinematic variables available - the and the invariant mass (M). 
While the spectra of dilepton is affected by the flow, the /ct integrated M spectra 
remains unchanged. This suggests that a careful selection of kr and M windows will 
be very useful to characterize the QGP and the hadronic phases. 

The interferometry of the dilepton pairs actually reflect correlation between two 
virtual photons, the analysis then concentrates on computing the Bose-Einstein cor- 
relation (BEG) function for two virtual photons which can be defined as, 

C (k k)-l I cos(Aa)]^ + [/ d'x uix, K) sin(Aa)]^ 



where ki is momentum of the individual photon, K = {ki + /c2)/2, Aa — a-i — 0:2, 

ttj = tMit cosh(|/j — v) ~ ^^iT cos{9 — ipi), MiT = \J kfrp + is the transverse mass 
and u:{x, K) is the source function related to the thermal emission rate of the virtual 
photons per unit four volume (see [3,4] for details). 

For the space time evolution of the system relativistic hydrodynamical model with 
cylindrical symmetry [5] and boost invariance along the longitudinal direction [6] 
has been used. For a system undergoing isentropic expansion, the initial tempera- 
ture (Tj) and proper thermalization time (t,) of the system may be constrained by 
the measured hadronic multiplicity, dN/dy ~ TfTi. Here we have taken Tj = 290 
MeV and Xj = 0.6 fm/c. The equation of state (EoS) which controls the rate of ex- 
pansion/cooling has been taken from the lattice QCD calculations [7]. The chemical 
(Tc/i=170 MeV) and kinetic (T/-o=120 MeV) freeze-out temperatures are fixed by the 
particle ratios and the slope of the kr spectra of hadrons [8] . With all these ingre- 
dients the correlation function C2 has been evaluated for different invariant mass 
windows ((M) = (Mi + M2)/2=0.3, 0.5, 0.7, 1.2, 1.6 and 2.5 GeV) as a function of 
qside and qout [3] which are related to transverse momentum of individual pair. Once 
the correlation function for the (time like) virtual photon is calculated, the source 
dimensions can be obtained by parameterizing it with the empirical (Gaussian) form: 

c,^i + ^eM-Rl<i)- (2) 



where the subscript i stand for out and side. In Eq. 2 A represents the degree of 
chaoticity of the source. Deviation of A from 1 will indicate the presence of non- 
chaotic sources. 

The Rside-i radius corresponding to qsi^e is closely related to the transverse size of 
the system and considerably affected by the collectivity and the Routi radius corre- 
sponding to qout measures both the transverse size and duration of particle emission 
[9,10]. The Rside shows non-monotonic dependence on (M) (Fig. la). It can be shown 
that Rside ~ ^/ {}- + Ecoiiectivel Ethermai)- With the trausvcrsc expansion the radial size 
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Fig. 1. The variation of Rside and Rout as a function of (M) (left panel). The right panel 



show Rout /Rside and Rdiff = \J Riut - ^lide ^ ^ function of (M). 

of the emission zone decreases with time as a rarefaction wave moves toward the 
center of the cyhndrical geometry. The high (M) regions are dominated by the early 
partonic phase [3] where the collective flow has not been developed fully [11] con- 
sequently the ratio of collective to thermal energy is small- hence a larger R^ide is 
obtained for large M. In contrast, the lepton pairs with M ~ nip are emitted from 
the late hadronic phase where the size of the emission zone is smaller due to larger 
collective flow giving rise to a smaller Rside- The ratio of collective to thermal energy 
for such cases is quite large, which is reflected as a dip in the variation of Rside with 
(M) around the p-mass region (Fig. la). Thus the variation of Rside with M can be 
used as an efficient tool to measure the collectivity in various phases of matter. The 
dip, in fact vanishes if the contributions from p and u is switched off [4]. We observe 
that by keeping the p and uj contributions and setting radial velocity, f,- = 0, the dip 
in -Rside vanishes, confirming the fact that the dip is caused by the large radial fiow 
of the hadronic matter [4]. The Rout probes both the transverse dimension as well 
as the duration of emission and unlike Rside-, Rout does not remain constant even in 
the absence of radial fiow. The large M regions are populated by lepton pairs from 
early partonic phase where the effect of flow is small and the duration of emission is 
also small - resulting in smaller values of Rout- For lepton pair from M ~ rup region 
the flow is large which could have resulted in a dip as in Rside in this M region. 
However, Rout probes the duration of emission too which is large for hadronic phase. 
The larger duration overwhelms the reduction of Rout due to flow in the hadronic 
phase resulting in a bump in Rout in this region of M (Fig. lb). 



J^nut ^side 



Figs. 1(c) and 1(d) show the variation of Rout I Rside and Rdiff = y -^^out 
with (M). These two quantities provides information on the duration of particle 
emission [10] for various domains of M. In one hand the high and low M domains 
are dominated by radiation from early QGP phase [4] where the duration of particle 
emission is small. On the other hand the lepton pairs for M around p mass dominantly 
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originate from the late hadronic phase where the duration of particle emission is large, 
Resulting in a non-monotonic variation of Rdiff, which indicates the presence two 
different phases during the evolution. 

In summary, we have evaluated the dilepton pair correlation functions relevant for 
Au+Au collisions at RHIC energy. The values of HBT radii extracted from the dilep- 
ton correlation functions show non-monotonic dependence on dilepton pair mass, 
reflecting the evolution of collective flow in the system. It appears that the non- 
monotonic variation of Rout and Rside with M originate from the presence of two 
phases during the evolution of the system, hence could be use as a signal of phase 
transition. 
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